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O Abstract 
CN 

Qh Five dimensional dilaton models are considered as possible holographic duals of the pure gauge 

(D 

C/^ QCD vacuum. In the framework of these models, the QCD trace anomaly equation is considered. 

00 

Each quantity appearing in that equation is computed by holographic means. Two exact solutions 
for different dilaton potentials corresponding to perturbative and non-perturbative /3-functions 
p I are studied. It is shown that in the perturbative case, where the ^-function is the QCD one 

at leading order, the resulting space is not asymptotically AdS. In the non-perturbative case, 

^ the model considered presents confinement of static quarks and leads to a non-vanishing gluon 

> 

condensate, although it does not correspond to an asymptotically free theory. Calculating the 

(N 

cr^ Nambu-Goto action corresponding to a small circular Wilson loop, leads to an expression for the 

gluon condensate. The validity of the trace anomaly equation is considered for both models. It 

O 

p^ holds for the perturbative model and it does not hold for the non-perturbative one. 
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I. INTRODUCTION 



The relation between large N gauge theories and string theory fT] together with the 
anti-de Sitter/conformal field theory (AdS/CFT) correspondence [2HH] have opened new 
insights into strongly interacting gauge theories. The application of these ideas to QCD has 
received significant attention since those breakthroughs. From the phenomenological point 
of view, the so called AdS/QCD approach has produced very interesting results in spite of 
the strong assumptions involved in its formulation [6]. It seems important to further proceed 
investigating these ideas and refining the current understanding of a possible QCD gravity 
dual. This endeavor has been followed in references [7]. The aim of the present paper is 
to explore the simplest non-perturbative features of QCD. This is done in the framework 
of a holographic description of the pure Yang-Mills (YM) QCD vacuum by means of 5- 
dimensional dilaton gravity models. 

At the basis of the AdS / CFT correspondence is the connection between scale transforma- 
tions in the boundary field theory and isometrics of the bulk gravitational theory. However, 
QCD is not a conformal field theory, as the scale symmetry is broken by the trace anomaly 
[8]. The trace anomaly equation describes the behavior of QCD under scale transformations. 
The question to be explored is how a holographic model can incorporate this behavior. 

The trace anomaly equation [8j states that, 

Ti = ^Tr{G,,G'^) (1) 

where denotes the trace of the QCD energy momentum tensor (latin indices for space- 
time), /3(A) is the QCD /3-function, X = N is the t'Hooft couphng, G,j is the QCD field 
strength tensor and the trace is taken in the fundamental representation of the SU(N) gauge 
group. In this respect it is important to note that holographic models can tell something 
about each of the three quantities involved in the trace anomaly equation, namely the 
vacuum expectation value (VEV) of the trace of the energy momentum tensor, the /3-function 
and the VEV of Tr(GyG''-'). 

According to the correspondence, evaluating the five dimensional action at a classical 
global solution gives information about the VEV of the trace of the energy momentum 
tensor. The /3-function can be obtained in terms of the solutions to the 5-dimensional 
equation of motion derived form the action in the bulk. Finally, there is a way of calculating 
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the VEV of the Wilson loop by means of minimizing the Nambu-Goto (NG) action for a 
loop lying in the boundary space. This is known to work in the strictly AdS case, i.e. 
for a conformal boundary field theory, and its generalization to non-conformal cases is still 
an open important problem. In turn, the VEV defined by G2 = ^^^^iGijG^^), known as 
the gluon condensate, can be determined from the coefficient of the area squared in the 
expansion of a small Wilson loop in powers of its area P-[TT]. 

The features and results of this work are summarized as follows, 

• Two exact solutions of 5-dimensional dilaton gravity for different dilaton potentials are 
considered. The first model, to be referred to as perturbative model, has a /3-function, 
which to leading order in the t'Hooft coupling is the same as the perturbative 1-loop 
QCD /3-function. The second model will be referred to as non-perturbative model 
(because its /3- function is non-analytic in A). This model, by choice of the parameter 
a in the model, can be made to correspond asymptotically to the soft wall model 
often used in non-dynamical models of holographic QCD. The model leads naturally 
to confinement in the sense of static quarks, and to a non-vanishing gluon condensate 
when tested with a Wilson loop. However, it does not lead to asymptotic freedom in 
the ultraviolet. 

• For the perturbative model the asymptotic behavior of the solutions in the ultraviolet 
is not AdS. In the language of the holographic renormalization group the difference 
with the AdS limit is produced by an irrelevant operator that fiows away from the 
AdS fixed point. In the non-perturbative model considered, the /3-function gives rise 
to an UV fixed point at finite A and the metric is asymptotically AdS. 

• Using the correspondence, the VEV of the energy momentum tensor is obtained by 
evaluating the 5-dimensional action in the corresponding exact solutions, regularizing 
by introducing an energy scale and subtracting. These subtractions are performed as 
proposed in [12], and employed in the holographic case in [l3j. In the perturbative case, 
taking into account Eqn. ([T]), it is argued that the same solution should be subtracted, 
leading to a vanishing VEV for the energy-momentum tensor. In the non-perturbative 
model, being asymptotically AdS, the AdS limit is subtracted. 

• In order to calculate the gluon condensate the VEV of a small circular Wilson loop is 
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considered. This is carried out using the corresponding NG action. For the perturba- 
tive model this procedure leads to a vanishing gluon condensate, while a non-vanishing 
result is obtained in the non-perturbative case. 

• The validity of Eqn.Q is considered for both models, and shown to hold in the pertur- 
bative one. In the non-perturbative model the dependence of the gluon condensate on 
the energy scale is not the one required by Eqn. ([T|. This is however not unexpected 
as this model does not give a consistent description of the QCD ultraviolet behavior. 

The paper is organized as follows. Section II presents the 5-dimensional dilaton-gravity 
model employed in what follows. Exact solutions of the dilaton model equations of motion 
and associated /3-functions corresponding to the perturbative and non-perturbative models 
are studied in section III. Section IV deals with the evaluation, regularization and subtraction 
of the gravitational action evaluated in the above mentioned exact solutions. Section V 
discusses the relevant asymptotics of the solutions of section III, and gives the explicit result 
for the subtracted gravitational action for those solutions. Section VI presents a study of the 
VEV of a small circular Wilson loop by means of the minimization of the NG action. Section 
VII addresses the issue of validity of the trace anomaly equation in the models considered. 
A final section VIII presents conclusions and outlook. 



II. DILATON MODEL 



The model considered is that of a self interacting scalar field immersed in a dynamical 
gravitational field in d + 1 dimensions (in the end the results are only valid at d = 4). The 
action of the model is given by [14], 

Sd+i = W(/ d^+\^g{-R+\g^''d,cl)dA-Vm-2 f d'xVhK], 

(2) 

where G^^^^^ is the Newton constant in (i+l-dimensions (of dimension (rf— 1)), g^^ the metric 
tensor field, R the scalar curvature, the dilaton field, and V{(f)) the dilaton potential. The 
last term is the Gibbons-Hawking term [15] where K is the second fundamental form. This 
term is included to make the Lagrangian depend only on the first derivatives of the metric. 
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The equations of motion derived from this action are, 

E^, - \ d^cj> d,ct> + ^-g^, {dcj>f - ^ ^^,1/(0) = 

d,{,/9 9'''d.<t^) + ^^^ = (3) 

where the Einstein tensor E^^, reads: Efii, — R^j,^, — \g^i,R, and (50)^ = g^''dn(j)d,y(j). Because 
here the focus is on the vacuum of the boundary field theory, only metrics and scalar fields 
having flat boundary space isometry invariancc are considered, thus only solutions for the 
metric and scalar field of the following general form are considered, 

ds^ = du^ + e^^^''^7]ijdx'dx\ (t> = (p{u), (4) 

where r]ij is a flat metric, and the coordinates employed here are known as domain wall 
coordinates. The boundary of the space is at m = ±cxo. The AdS metric corresponds to 
taking AAds{u) — u, where the coordinate u is measured in units of the AdS radius L. For 
this particular choice of flelds which only depend on u, the equations of motion are given 

by, 

A" + dA''-^ = 
o — 1 

dA'^ - — — - - = 

2{d-l) d-1 

0" + d^>' + ^ = O, (5) 
d(p 

where the prime denotes derivation with respect to u. Introducing a superpotential 
according to: 

A'{u) = W{(f)) (6) 
the choice ^ — 2{1 — d) < reduces the three equations in Eqn. (5) to the single equation: 

Since the intended realistic application to QCD is at d = 4, throughout ^ = — 6 could be 
replaced. 
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III. /3-FUNCTIONS IN DILATON MODELS 



In the AdS/CFT correspondence the identification is made of the YM couphng with the 
dilaton profile according to [5]: 

The energy scale /i (measured in units of a scale where L is the length unit mentioned 
earlier) of the boundary theory is identified with the scale factor e"^'-"-' in domain wall coor- 
dinates: jji = e^^^\ These identifications give the /3-function in the dilaton model [?]: 

In the rest of this section two different and exactly soluble dilaton models are considered. 
These models are obtained according to the following scheme: a dilaton profile is given, 
where by expressing (p'{u) in terms of (p{u) and employing Eqn. ^ the superpotential W{(j)) 



is obtained, followed by integrating Eqn. (|6]) to obtain A{u), and finally from Eqn. (10) the 
/3-function is obtained. The potential V^(0) is determined from Eqn. (|5]). 

The two models considered are extreme cases. One model corresponds at leading order 
in the gauge coupling to the perturbative QCD /3-function, while the other one corresponds 
to a non-perturbative /3-function, i.e., which is non-analytic at small coupling and which 
leads to an UV fix point. These models are qualitatively different as the next sections show. 
The precise choice of dilaton profiles is made so as to be able to perform all the calculations 
analytically. 

A. Perturbative /3-function 

The following dilaton profile is considered, 

^(u) = -Uogi{auf + K') (11) 

Note that this choice means that A < N/k. Therefore, k should be a quantity order N. 
Using the procedure just described leads to: 

A{u) = Ao + Aiu-— log {{a uf + k^) + — — arctan — ) , (12) 



6 



where for convenience the integration constants can be chosen in such a way that the leading 
asymptotic behavior be AdS, namely Aq = j^, and Ai = 1 — f^. Then, asymptotically 
A{u) ~ M — I log(a'u) + 0{l/u^). The resulting /3 function reads: 



which to leading order in A becomes 



/3(A) = ^ ^, (13) 



m = -^ + o{x'). (14) 

The choice a = reproduces the leading order term of the QCD /5-function (see Fig. 1) 



B. Non-perturbative /3-function 

A /3-function with non-perturbative behavior, i.e. non-analytic in the coupling A, is 
obtained from the following dilaton profile, 

0(M) = Ce-"" (15) 

where a > 0. In this case, 

A{u)=u+^e-'-, (16) 

giving an asymptotically AdS metric. 

The resulting /3-function is then given by. 



« A log 4 , , 

which is positive in the interval < A < A^, leading to an UV fixed point at A = A^ (see 
Fig. 1). Thus, this theory is not asymptotically free, and therefore is not related to a pure 
YM theory. The sign of the constant C determines two phases of the theory: for C < the 
theory becomes free in the infrared, while for C > it becomes strongly coupled. Indeed 
this latter case describes a confining theory in the IR. In order to see this it is convenient to 
express the above result in the conformal coordinate z, where asymptotically u = — log(z), 
and therefore A(z) = — \og{z) + ^-2^°- This A matches the Giirsoy-Kiritsis [7J criterion for 
confinement Pl] The negative sign of the coefficient multiplying the z'^°' term is crucial in 
two respects: it is necessary for the confinement criterion [7j to be fulfilled and second, the 




X 



FIG. 1: P functions of the models considered: perturbative (dashed), non-perturbative (black), 
perturbative QCD (gray). (5 = corresponds to the AdS case. 

behavior of the factor e^^^^ for 2; — > oo is such that, lim2_>.oo e^*^^) = 0, which as shown in 
the next section, makes the use of a infrared cut-off unnecessary in the evaluation of the 
5-dimensional action for this solution. 



IV. THE TRACE OF THE ENERGY-MOMENTUM TENSOR 

According to the AdS / CFT conjecture, taking the metric as the source field of the energy- 
momentum tensor of the boundary field theory, the VEV of the trace of the energy momen- 
tum tensor is evaluated by simply evaluating the action Eqn. (2) for the classical solutions 
of the previous section. 

Taking the trace in the first Eqn. (3) gives. 



R 



jd+l) 
(l-d) 



v{<p) + -idci>r, 



(18) 



and the action for the classical solutions becomes: 




(19) 



and using the first Eqn. (5), 




(20) 
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Noting that J^e'^^(") = d e'^^(")(A" + dA'^) and ^ = e-^^^") leads to 

Sbui, = r du ie^^W 

16 



1671 G^i"-'^ d du^ 



^^dA{u) 

du 



(21) 

boundary 



where Vm^ denotes the volume of the boundary d- dimensional space. On the other hand, 
the classical Gibbons-Hawking boundary action is given by: 

Sgh = [ d'xVh, (22) 

where h is the induced metric in the boundary M^, namely y/h = e'^"^^'^\ and denotes 
a unit vector field orthogonal to the boundary of M^+i- In domain wall coordinates this 
vector field is simply -i- = and therefore: 



dn du 

d 



^GH — TJZT^^Mj^ 



du 



^dA(u) 



(23) 

boundaTy 



'N 

which is just —d times the bulk action. For both exact solutions considered in the previous 
section there is no contribution from the infrared boundary. On the other hand, the ultravi- 
olet boundary u ^ oo gives for both cases divergent contributions, as it happens in general 
for any holographic model. As proposed in these contributions can be regularized by 
evaluating at a finite value Uq. This leads finally to, 

S = — ^1^M,(1 - rf)e'^^(-)A'(t.o). (24) 

It is important to note that for a boundary theory that is not quantum conformal invariant, 
as for example QCD, the regulator uq has a physical meaning. Indeed, as mentioned in the 
previous section, the energy scale at which the boundary theory is observed is related to mq, 
the boundary value of the domain wall coordinate u. 

As shown in p!2j and applied to holographic models in [13], a well defined action can be 
obtained by subtracting from the regulated action an action corresponding to some back- 
ground metric having the same asymptotic limit. That is, 

Ssub. = S-S'''y'^''-, (25) 

where S"''^^^^' denotes the action evaluated in a solution having the same asymptotic behavior 
as the classical one. The subtracted energy-momentum tensor is obtained recalling that, 



according to the correspondence, 

S = [ (fx VhhijT\ (26) 

J Ma 

leading to Tl{sub.) = , where A{uo) denotes the common asymptotic exponent. 

The choice of this background metric for the solutions considered in section III is discussed 
in the next section. 



V. THE UV QCD FIXED POINT 

The perturbative model in subsection IIIA presents features the understanding of which 
leads to new insights. These are the following: 

• The model leads to a /3-function that coincides at leading order with the perturbative 
QCD /3-function. 



• The model is not asymptotically AdS. As Eqn. (12) shows, the deviation of A{u) from 
the AdS limit becomes ^ logu. 

• As shown in the previous section, the action should be subtracted with the action 
evaluated in a background metric having the same asymptotic behavior as the one 
to be subtracted. Thus, it is not sufficient to perform a substruction with the AdS 
metric. 

• In the language of the holographic renormalization group [T7], this correction corre- 
sponds to an irrelevant operator, that flows away from the AdS fixed point [5J. This 
can be seen from the fact that the dilaton field behaves as — log u at the UV boundary. 

• Eqn. ([T]) implies that for QCD the trace of the energy-momentum tensor should vanish 
in the UV. This can be independently seen in two ways. As shown in section VI for 
this model, the VEV of the Wilson loop, calculated via the NG action, does not have 
terms which are powers of its area, and therefore the gluon condensate G2 must vanish. 
The other way is simply to recall that in perturbative QCD the log of the VEV of the 
Wilson loop follows a perimeter law. 

All these points indicate that the UV fixed point of QCD does not correspond to AdS. 
It corresponds to another solution that is well approximated by the one in subsection IIIA 
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in the UV, i.e. for large m, and therefore the action evaluated in the same solution or one 
asymptotically equivalent must be subtracted, leading to a vanishing trace of the subtracted 
energy-momentum tensor. 

In the non-perturbative model the space is asymptotically AdS, and the subtracted action 
becomes: 



Sfl = - S^^' = fe#-^"-(l - ^e---) - 1 ) , (27) 



leading to, 



^2 _2 



THs,A, NP) = (l - e-". - e-« ) (28) 



VI. WILSON LOOPS 



The VEV of the operator G2 (gluon condensate) appearing in the trace anomaly is ac- 
cessible through the power like behavior of small Wilson loops as a function of their size. In 
pure YM theory the expansion of a small smooth Wilson loop (e.g., square or circular) is 
expected to have the form given by [HHUl [IS] : 

n) "T2iv^^'' + --- ^^^^ 

where i is the length of the loop, s is its area, and Z = /3i(A)//3(A) where (3i is the one loop 
/3-function. It is argued in pure YM that the terms order s vanish as these would require a 
gauge invariant dimension two condensate. 

The connection between Wilson loops of the boundary conformal gauge theory and min- 
imal surfaces was made in references [TSl [12]. According to it, in a CFT such as A/" = 4 
SUSY YM, in the large N limit and large 'tHooft coupling the VEV of the Wilson loop is 
determined by the minimal area surface in the d + 1 AdS space subtended by the path of 
the loop r. Specifically: 

W{T) = ^ Tr PExp(- ^ Aidx') 
{W{V)) oc e-^^ (30) 

where the minimal area S'r is given by the NG action of a string whose ends run along the 
loop. Since for a loop located at the boundary S'r diverges, it has to be regulated, and thus 
the proportionality factor above. 
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The extension of this identification to non-conformal YM theory is still an open problem, 
in particular because in that C3iSG, clS discussed earlier, the theory cannot be obtained via a 
relevant deformation of a CFT [5J . This problem is closely related to the problem of finding 
the non-critical string action for QCD [20j. An extension of the correspondence for Wilson 
loops to the non-conformal case has been proposed [7j, in which the NG action is the one 
corresponding to the string-frame metric, namely: As{z) = A{z) + 0(2)/v^ in c? = 4 
dimensions. 

For the present purpose a circular Wilson loop of radius a is considered, for which the 
NG action turns out to be: 

Sng = f'dpp e2^^('^-('')) v^TWG^, (31) 

2tt a' Jq 

where r = a p is the radial coordinate of the disk, and z = a u is the bulk coordinate in 
conformal coordinates. The equation of motion is: 

pcu" + {I + uj'^){uj' - 2a p A'siatu)) =0, (32) 



where the solution needed satisfies uj{1) = 0. In AdS limit it is ui{p) = a/1 — p^, a half 
sphere. 

The UV divergencies of the NG action result from the contributions to the integral for 
p — > 1. Noticing that w'(p) diverges as p — t- 1, one obtains: 

dS^G _ ae'^sM 
ozq 271 a 

where zq can be interpreted as the location of the loop in the bulk coordinate z (provided 
Zq << a). In dilaton models one readily obtains: 

dSMG ^ ae'^s(zo)-Ao 

dAo 2TTa'W{Ao)' ^ ^ 

where Aq = A{zq), which asymptotically for the models discussed A{zq) — ^ — log^o- If the 
/3-function is given as input to the model, the superpotential and A{z) are given by 



W{X, = exp (i / ^dX 



MX, - (35) 

One readily checks the AdS case where W = const and As = A, giving ^^§j^ = 2^ 
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The perturbative model asymptotically gives ^(A) = —^X^, — — log(Q;^) and 

= 1 — where, without loss of generality, the constant of integration required for 
Vr(y4) has been chosen to be Wq — 1. This leads to: 



- exp (^Aq - log(aAo; 



(36) 



dAo 271 a' 

This shows that, as one would expect from the fact that the metric is not asymptotically 
AdS, the UV divergence of the action is modified with respect to the AdS case by the second 
term in the exponent. 

The non-perturbative model is asymptotically AdS and thus the expectation is that the 
UV divergence coincides with the AdS case. If the coefficient a > 1 this is indeed the case 
as it is easily shown using Eqns. (15) to (17) for a > 1, which leads to: 

For a = 1 a constant term remains, which corresponds to a term linear in Aq in the UV 
divergence of Sfqc or equivalently logarithmic in zq. 

The UV divergencies stem from the fact that As diverges at the boundary. Therefore, 
they must naturally be only proportional to the perimeter of the loop, i.e. proportional to 
a. For this reason, the contributions of higher powers of a, which are of interest here, are 
independent of the regularization of Sjsig and unambiguous. 

The central point of the discussion is the sufficient conditions for the presence of higher 
power terms in a in S^q- The simplest case is when the metric is asymptotically AdS 
and the UV divergence of S^g corresponds as well to the AdS case. For small a, As{z) = 
^Adsiz) + SA{z), and expanding in SA leads to: 

S^G^^, f%{l + 2 M(ar?)) + 0{5A'), (38) 



where rj = y^l — p^, and evaluation in the AdS limit solution has been performed. The first 
order approximation is adequate near the boundary rj ^ only if the UV divergencies are 
strictly AdS. On the other hand, the dependencies of S^g in powers of a beyond the first 
power (perimeter terms) will stem primarily from the interior of the integration domain, 
where the approximation is expected to work. Thus, a sufficient condition for such power 
corrections is that 5 A contains terms which have power dependency in the argument. The 
contributions 0{5A'^) in Eqn. (38) are in general difficult to evaluate as they involve the 
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corrections to the solution of the equation of motion (32) [16]. The arguments made here 
apply in particular to the non-perturbative model when a > 1. 

When the metric is not asymptotically AdS, as is the case of the perturbative model, a 
more accurate evaluation is necessary. For a sufficiently small the entire surface will lye near 
the boundary m — )■ cxd, and k can be set to zero, thus </>(«) = — log(au), A{u) = u—^ log(au). 
Setting u = — log 2; and evaluating Sng with the asymptotic AdS solution 2(77) = ar] leads 
to: 

Sng = ^ ^^^P^-^^ + 7f ) log(-«log(ar/))). (39) 
It is readily checked that this has the UV divergence obtained earlier in Eqn. (36). Evidently 
the dependence of Sng in cl is logarithmic, and therefore according to the evaluation of the 
Wilson loop G2 = in the perturbative model. A similar conclusion results if /3(A) is in 
general analytic in A. Therefore, in the present framework, this indicates that in order to 
obtain a non-vanishing gluon condensate, the /3 function should include non-analytic terms 
in A. 

As an illustration of the latter, where power corrections are obtained at small coupling 
as consequence of non-perturbative terms in the /3-function, consider the asymptotically 
free theory with /3(A) = — 6oA^(l + c exp(— ^)), which is found in certain SUSY gauge 
theories [21] as the result of instanton contributions. Considering the non-perturbative 
piece as small (or expanding in c), asymptotically H^(y4) = e~«^(l — ^^^^^^2 6""^""^), 0(^) = 
- logboA + ^e-"^«^, leading to: 

dSNG _ OSng'^''-. 2c^^ 1 + aol 

which as expected coincides asymptotically with the perturbative model. The evaluation of 
the finite pieces gives power terms in a. Asymptotically, to first order in c: 

As = Al'^'iz) + ^ (^-^ exp(-2a 60 Apert(^)) + exp(-2a 60 A^,,,{z))^ , (41) 

where pert, indicates the case with c = discussed earlier. Using Eqn. (38) leads to: 



Sng = ^^G* + ZTTT^ / S ( r\ . + iavr''] 



d7] (ar7)"*o 
ira'abo Jo r]^ y/Slogiar]) 

5^P°J- = \ -a'"'\ (42) 

a'abo [abo - 1) ^ ^ 
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obtained after replacing Apert ~ AAds in the evaluation. Note that the power correction 
in this case did not stem from the contribution to As by the dilaton, but rather from the 
correction order c to the metric A itself. This model gives a non-vanishing G2 if a&o = 2. 

The non-perturbative model is now analyzed for a > 1, where (f){u) = Ce^"", A{u) = 
u + ^e^"", and asymptotically u = — log 2;. Applying Eqn. (39) leads to: 

Sng = SMG{AdS) + C% {^{avT + 0{ar^f'^), (43) 

Jo 1] a/3 

where the term oc {at])"' stems from the contribution to As by the dilaton. Clearly, if a = 4 
the model gives a non- vanishing G2, namely G2 = a = 1 it reproduces the 

additional logarithmic contribution in zq to the UV divergence in Eqn. (38). For a = 2 the 
model is similar to the one analyzed in [16]. In that case, to obtain the power correction 
it is necessary to calculate to second order in the perturbation to the action, and therefore 
corrections to the solutions are to be calculated. As mentioned earlier, in QCD the power 
series in the area s of the Wilson loop start at ~ a^; for a = 2 there is however a 
non- vanishing term order 0? |16j . 



VII. THE TRACE ANOMALY TEST 



For the perturbative case the trace anomaly equation is clearly fulfilled. Indeed the 
subtraction to the 5-dimensional action in section V was performed in order to match, 
through Eqn. (1), the vanishing of G2 determined in the previous section for this model. 
On the other hand, for the non-perturbative case, it is shown below that it is not possible 
to match both sides of Eqn. (1). 



A. The trace anomaly equation for the non-perturbative case 



Equations (17) and (28) for a > 1 and ci = 4 lead asymptotically for — t- 00 to: 



/3(A) a<p aCe-""« aGz^ 



A (l + f|02) (l+^^e-2-o) (l + ^^g^2a) 

Tl{suh.,NP) = --^(l-e^-^""" + ^e---) 



(a + 2)C^ 



2a 



2,2tiG^ 



(5) ^0 ' 



(44) 



(45) 
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where the asymptotic relation between domain wall and conformal coordinates zq = e~"° has 
been employed. If the trace anomaly equation Eqn. ([T]) would be fuUfiUed, replacing Eqns. 
(44) and (45) into Eqn. ([T]) would imply for the gluon condensate to vanish asymptotically 



as: 

a + 2 



G,{zo) = {Cz- + C'ztr + ■■■)• (46) 

B. Wilson loop calculation of G2 

The computation of G2 using the Wilson loop calculations of the previous section involves 
a different choice of boundary conditions than the one employed in this section. This is 
because the Wilson loop should be situated at a finite value of the coordinate orthogonal 
to the boundary, corresponding to the finite value chosen in evaluating the 5-dimensional 
action used to evaluate the trace of the energy-momentum tensor. The boundary condition 
to be employed is, 

z{a) = zq. (47) 
For the pure AdS case, a solution of the area minimization equation satisfying this boundary 



condition is given by, z{r) = a/o^ — + Zq, which simply corresponds to a circle of radius 



R = a'^ + Zq that is the radius required to match the boundary condition (47). For the 
non-perturbative model the effect of the above mentioned change in boundary conditions 
is well aproximated by replacing the radius a by the effective one corresponding to the 
AdS solution, i.e. R = \/a?~-\-~z^. Making that replacement in Eqn. (43) shows that the 
coefficient of a" has a contribution, coming from the term proportional to i?" [25j, which 
does not vanish for zq = 0. In particular, the simplest case where a = 4 gives a putative 
G2 7^ 0. This is however in contradiction with the dependence in Eqn. (46), which comes 
from assuming the validity of ([T]). Therefore, the trace anomaly equation is not fulfilled in 
this model, and this is so in general for a > 1. 



VIII. CONCLUSIONS AND OUTLOOK 



In this work the validity of the trace anomaly equation has been studied in the holographic 
framework. This was done by considering holographic evaluations of the VEV of the trace of 
the energy-momentum tensor, the /3-function and the gluon condensate 6*2- The /3-function 
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is directly related to the definition of the particular model under consideration. The VEV 
of the trace of the energy momentum tensor was evaluated according to the holographic 
correspondence, by evaluating the d + 1 dimensional classical action of the dilaton model 
on the corresponding classical solution. The gluon condensate can be obtained in a YM 
theory from the VEV of the Wilson loop, which was here evaluated for the models studied 
by means of a NG action. 

Two models were analyzed, which can be exactly solved and which have different qual- 
itative characteristics. In the perturbative model, where G2 = 0, consistency is fulfilled as 
the evaluation of the classical action can be appropriately subtracted to give a vanishing 
trace for the energy momentum tensor. If indeed G2 = in QCD, this may already be a 
somewhat realistic model. On the other hand, the non-perturbative model shows an incon- 
sistency for the trace anomaly equation. This is manifested by the fact that G2 has different 
behavior in the scale zq in the two evaluations. Indeed, the evaluation of the action and 
Eqn. ([1]) give G2 oc e~""° = Zq, while from the Wilson loop evaluation G2 is non-vanishing 
in the limit Zq — )■ 0. This inconsistency seems reasonable since the non-perturbative model 
fails to correctly describe the UV properties of QCD, being asymptotically AdS and not 
asymptotically free. 

Various interesting conclusions can be drawn from these results. They indicate that, 
although a holographic model of the pure gauge QCD vacuum based on the AdS space is 
not feasible, they do not preclude a gravitational dual based on a dynamical 5-dimensional 
Einstein gravitational theory. They also show that QCD Ward identities, as for example the 
trace anomaly equation, strongly restrict the possibilities. It is reasonable to expect that 
QCD symmetry restrictions can in principle lead to a more precise version of its putative 
gravitational dual. Such a dual should lead to a boundary theory having all the following 
properties: asymptotic freedom in the UV, confinement in the IR, (possibly) a non-vanishing 
gluon condensate, and consistency with the trace anomaly equation. As the examples con- 
sidered have shown, it is not at all obvious how to obtain a consistent model with these 
properties. Work in this direction is in progress and will be reported in due course. 

Among important fundamental non-perturbative effects in QCD, the existence of a non- 
vanishing gluon condensate was early on identified [22] . It has important manifestations in 
hadron phenomenology [2S], and there are indications of its non- vanishing from lattice 
QCD pm ITT] . Due to its importance, its further understanding in the framework of holo- 
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graphic models of QCD is going to play a key role in the development of such models, as it 
has been shown in this work. 



IX. APPENDIX 

This appendix presents an exphcit calculation of the NG action in a case where the As{z) 
has deviations from the AdS limit which are integer powers of z, namely, 

As{z) = - log z + Y.c^nz''. (48) 

n 

The equation of motion Eqn. (32) is solved using an asymptotic series: 

n 

coiv) = 7j{l + J2J2Cne r^-logS), (49) 
n e=o 

where the coefficients Cneica, a) are obtained in a systematic fashion. 

The evaluation presented here can be applied to the non-perturbative model discussed in 
the text. A straightforward but lengthy evaluation gives: 

= 2W [---^<^^i^ogizo/a) + -aa^ 

+ ^3.32435 + y ^2 
+ (3.12395 + 5.03896 eta ai + 2 ag) (50) 
+ 1^12.4174 al + 19.5861 as + 2.09778 "i + 6.4849 "2 + y "4^ + C'(a^) 

For instance, in a "soft wall" model where only q;2 7^ one obtains: 

^soft wan(^) ^ ^ ^^^2(_5^ + rj' + ...) + aWi-^-^v + ^-^v' + ■■■) + ■■ ■) (51) 

and the resulting NG action becomes: 

C-soft wall 1 I ^ ^2 , 134821 2 4 , \ /rn\ 

^ w + y + ^0790""^'^ + • • • ^ ^^^^ 

For the non-perturbative model with a = 4, one keeps only the term with 7^ 0, and the 
NG action becomes: 

^- = 2i?(^ + ¥"'»' + ^"'^0- '''' 
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